Terahertz conductivity spectroscopy of Co-doped BaFe 2 As 2 Thin Film 
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We investigated the complex conductivity spectrum of a Co-doped BaFe2As2 epitaxial thin film in 
the THz region. In the normal state, the complex conductivity shows a Drude-type frequency depen- 
dence, while in the superconducting state, the frequency dependence of the complex conductivity 
changes to that of a typical superconducting materials. We estimated the magnetic penetration 
depth at absolute zero to be 710 nm and the superconducting gap energy to be 2.8 meV, which 
is considered to be the superconducting gap opened at the electron-type Fermi surface near the M 
point. We succeeded in obtaining the low-energy elementary excitation of a Fe-based superconductor 
using the electromagnetic method without invoking the Kramers-Kronig transformation. 



Fe-based superconductors have stimulated significant 
interest and competition for the synthesis of a novel type 
of superconductor since the discovery of superconductiv- 
ity in F-doped LaFeAsO in 2008 -. Currently, Fe-based 
superconductors are classified into five families based on 
the structure. Among them, i?FeAsO (R : rare earths) 
and j4eFe2As2 (Ae = Ca, Sr, Ba, Eu, K) systems have 
been widely investigated^. The highest superconducting 
transition temperature, T c ~ 55 K, is found in F-doped 
i?FeAsO£. Yet, the AeFe2As2 system is superior for thin 
film application, due to its relatively high T c and the 
existence of only one kind of anion, which makes the fab- 
rication of the high-quality thin film easier than in the 
i?FeAsO system. 

As for the AeFe2As2 system, the undoped BaFe2As2 
possesses metallic conductivity^, and Co-doping of the 
Fe-site leads to electron-doped superconductivity at 22 
K£. With decreasing temperature, structural and an- 
tiferromagnetic phase transitions occur at 134-140 
that are accompanied by an anomaly in the DC resis- 
tivity. These phase transitions are observed by the DC 
resistivity, up to 6 % Co-doping in the electronic phase 
diagram&&£. Superconductivity occurs at 2.5 - 18 % Co- 
doping, which indicates that there is a subtle relationship 
between the antiferromagnetic phase and the supercon- 
ducting phase in the underdoped BaFc2As2- The Fermi 
surface of Fe-based superconductors is much more com- 
plicated than that of high-T c cuprate superconductors 
because of their multiband nature. A band calculation 
study clarified that there are hole-type and electron-type 
Fermi surfaces around the Y and M points in a reciprocal 
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space, respectively^. The superconducting gap opens at 
both Fermi surfaces in the superconducting state, which 
causes the novel symmetry of the superconducting gap, 
s±— . By the thermal conductivity^ 2 -, the microwave 
conductivity-^ 3 -, and the ARPESi^, the existence of the 
nodeless gap in the superconducting state is suggested, 
which is consistent with the s± symmetry. The Fermi- 
surface-dependent superconducting gap energies, which 
are identified in a Ba(Feo.925Coo.o75)2As2 single crystal 
(T c = 25.5 K) by ARPESii, are 6.7 meV at the 'T" 
surface and 4.5 meV at the "M" surface, respectively. 

For the fabrication of AeFe2As2 thin films, we suc- 
ceeded in obtaining Co-doped SrFe2As2 epitaxial thin 
films 15 . Because the films are quite sensitive to water 
vapor in ambient air—, it is not easy to measure their 
superconducting properties. However, we recently found 
that Co-doped BaFe2As2 thin films are much more stable 
in air than the Sr-based thin films 17 , which increases the 
possibility of studying their physical properties, including 
the superconducting gap energy, A sc . The energy scale 
of A sc is close to the terahertz (THz) frequency range, 
thus we investigated the THz conductivity of Co-doped 
BaFe2As2 thin film in this paper. 

Co-doped BaFe2As2 epitaxial thin film was grown on 
a MgO(OOl) substrate at 730 °C by a pulsed laser depo- 
sition technique using a BaFe1.sCoo.2As2 disk target—. 
The overall dimension of the film was 5x10 mm 2 , and 
the film thickness was 110 nm, which was measured us- 
ing a surface profiler. Figure 1(a) shows the temperature 
dependence of the DC resistivity. There is no anomaly 
indicating the antiferromagnetic phase. From the en- 
larged plot near T c (the inset of Fig. 1(a)), the onset 
temperature of the superconducting transition, T° nsct , 
is estimated as 19.9 K. The midpoint of the transition 
is 18.7 K, and the transition width is approximately 
2.9 K. The X-ray diffraction measurements indicate that 



2 



0.5 
^ 0.4 

a 

<-> 0.3 

a 

^ o.i 
o.o 



Co-doped BaFe 2 As, film 
7-™= 19.9 k 



a 

So.] 



(a) 




10 15 20 25 30 

r(K) 



50 100 150 200 250 300 

r(K) 




-1 1 

Acj(deg) 



A/2 



780nm Laser 



| pump beam | 

chopper 



2 kHz 



(c) 



p-type 
InAs 



| probe beam | 



DryN 2 



BS 



sample 

5-300 K 



time delay 



I i PC switch 

^C^^j THzwave| ^P ::= ^ K ^g^ 



pre- Amp *10 



Lock-in-Amp 



multimeter 



He 
compressor 



FIG. 1: (color online), (a) The temperature dependence for 
the DC resistivity of Co-doped BaFe2As2 thin film. The inset 
figure shows the enlarged plot of the region near T c . The 
onset temperature of the superconducting transition, T° nact , 
was determined from the crosspoint of the dotted additional 
lines, (b) The rocking curve of the (004) reflection, (c) The 
schematic figure of the transmitted-type THz time-domain 
spectroscopy. Details are described in the text. 



the epitaxial relationship between the film and the sub- 
strate is BaFe 2 As 2 (001) || MgO(OOl) for out-of-plane, 
and BaFe 2 As 2 [100] || MgO[100] for in-plane. Figure 1(b) 
is the rocking curve for the (004) reflection of the film, 
while the full width at half maximum is about 1 deg. 

The frequency-dependent (0.2-1.5 THz) complex con- 
ductivity was measured by the transmitted-type THz 
time-domain spectroscopy (THz-TDS) 18 , using the mea- 
surement system shown in Fig. 1(c). For the emitter 
and the detector of the THz-wave, we used a Zn-doped 
InAs crystal and a dipole-type photoconductive switch, 
respectively. By measuring the photocurrent of the THz 
detector as a function of the optical path length of the 
probe beam (blue line in Fig. 1(c)), we can obtain the 
time-domain waveform of the THz pulse. Two time- 
domain waveforms of the transmitted THz-wave through 
the sample (Co-doped BaFe 2 As 2 film on a MgO sub- 
strate) and the reference (MgO substrate only) were mea- 
sured. In the cryostat containing quartz optical windows, 
we set the sample holder such that the direction of the 
electric field was parallel to the FeAs plane. The work- 
ing surface area of the film was approximately lxl mm 2 , 
where the THz pulse was irradiated. The measurement 
system was continuously purged with dry nitrogen to 
eliminate the effect of the absorption spectrum of water. 

We obtained the frequency-dependent transmittance, 



T(uj) (uj is the angular frequency), and the phase shift, 
A(f>(tj), by dividing the Fourier-transformed signal of the 
sample by that of the reference. The complex refractive 
index of the film, n = n — in, was obtained by fitting the 
data to the analytical formulaic. 
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where n su b is the complex refractive index of the MgO 
substrate, c is the speed of light, and d is the thickness 
of the film, respectively. Then, the complex dielectric 
function, e = e\ — ie 2 = Son 2 , where £q is the per- 
mittivity of a vacuum, and the complex conductivity, 
a = (Ti+i(T 2 = iu>(s — £o)j were calculated. We separately 
measured n S ub for the MgO substrate as n su t> ~ 3.18 and 
K su b ~ 0.006, which are consistent with the previously re- 
ported values 2 ^. Therefore, h(ui) can be obtained with- 
out using the Kramers-Kronig transformation. For the 
THz-TDS method, the reliable amplitude and the phase 
can be obtained in a single experiment, without perform- 
ing numerrous experiments over a wide frequency interval 
as an usual optical measurement. 

In the transmitted-type THz-TDS, the difference in 
the thickness of the substrate, Sd, causes a serious er- 
ror in A(f>. To eliminate the error in the phase shift, 
we corrected the obtained phase shift in the following 
manner. Because n su b 3> K S ub, Sd has almost no in- 
fluence on the transmittance. If we assume Sd = 20 
lim, the error in the transmittance is roughly estimated 
by 1 - exp[-2cjK sub <5<i/c] = 0.005 at 1 THz. There- 
fore, the error from the experimentally obtained trans- 
mittance can be ignored, and only the experimentally 
obtained phase shift, A0 exp , requires calibration. In ad- 
dition, we assumed that n = k at the highest tempera- 
ture (250 K). This indicates that <r 2 = ujeo(n 2 — n 2 ) = 0. 
This corresponds to the low-frequency response of the 
Drude metal (w < r _1 <C lo p ) and is widely used for 
the discussion of the microwave conductivity of metal- 
lic materials^. Therefore, from the experimentally ob- 
tained T{uj) and n — k = 0, we calculated (n + k, A(f> n — K ) 
from eq. (1) at each frequency. Then, we extracted Sd as 
(A(j) exp - A(j) n=K )c/ ujn suh , and subtracted A^ oxp -A0 n=re 
from A(f> exp over the entire temperature region. Sd was 
then determined to be 23.51 ± 6.6 /im, which is consis- 
tent with the value roughly measured by the micrometer. 
Therefore, we believe that our treatment is appropriate. 
After this calibration, we calculated the a(u>) of the film. 

Figure 2(a) shows the temperature dependence of the 
amplitude of transmitted THz pulses through the Co- 
doped BaFe 2 As 2 film at T = 5.1-250 K. The inset of Fig. 
2(a) shows the time-domain waveforms of the THz pulse 
through the reference and the sample. In the normal 
state, the amplitude of the transmitted THz pulse de- 
creases with temperature, which is based on the reduc- 
tion of the DC resistivity. In the superconducting state, 
the amplitude of the transmitted THz pulse is suppressed 
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FIG. 2: (color online), (a) The temperature dependence of the amplitude of the transmitted THz waveform through the sample. 
The inset shows the time-domain waveforms of the THz pulse through the reference (thickest solid line) and through the sample 
(second thickest line: 250 K, thinnest line: 5.1 K). (b) The frequency dependence of the complex conductivity. The lower and 
the upper panel represent the real and the imaginary parts, respectively. The solid line in the bottom panel represents the 
DC conductivity (3xl0 3 fi _1 cm _1 ) at 200 K, while the solid line in the upper panel indicates the lj^ 1 fitting curve, (c) An 
enlarged plot of (72(a)) at T > T c . The solid lines are linear fitting curves. The inset is the temperature dependence of the 
lifetime of the carrier, which is estimated from the slope of (72(01). 



to a great extent. These types of temperature depen- 
dence are found in many known superconductors^. Fig- 
ure 2(b) shows ct{uj) (lower panel is the real part, while 
upper panel is the imaginary part). Under conditions 
where T > T c , o\ ~ ctdc (°"dc is the DC conductivity, 
which is represented by the solid line (at T = 200 K) in 
the bottom panel of Fig. 2(b)), and is almost indepen- 
dent of frequency. In addition, (72 is almost equivalent 
to zero and slightly increases with frequency. There- 
fore, g(uj) resembles that of Drude-type metal. In Fe- 
based superconductors, there are two kinds of carriers 
(holes at the Fermi surface near the T point and elec- 
trons at the Fermi surface near the M point). Therefore, 
the experimentally obtained conductivity is the summa- 
tion of the conductivity of both the electrons and the 
holes. However, for the electron-doped system, the vol- 
ume of the hole pocket decreases with increasing dop- 
ing level, in contrast to that of the electron pocket. For 
BaFe1.gCoo.2As2, the number of the electrons is approx- 
imately twice as large as that of the holes^ 2 -. Therefore, 
we estimate that the life-time of the electron is reflected 
in the measured value 02(0-'), as r = 02(0;) / 'ctdc^ which 
is shown in the inset of Fig. 2(c). The obtained plasma 
frequency, uj p — v^dcT^oT", is about 0.56 eV at T = 200 
K, and this frequency does not almost depend on the 
temperature from 50 K to 200 K. 

As the temperature decreases to a value below T c , the 
ct{lo) factor of a typical superconducting material can 
be observed. o\ is suppressed over the entire frequency 
region and 02 is significantly enhanced compared with 
that in the normal state. We note that there remains 
some finite conductivity (~ 4 x 10 3 J7^ 1 cm _1 ) in the 
superconducting state, which is qualitatively similar to 
the result of an optical measurement. We have some 
possibilities for the origin of the finite conductivity in 
the superconducting state. (1) the residual conductiv- 
ity by uncondensed carriers not because of the presence 



of nodes but because of the dirty nature of the film we 
used, (2) the residual conductivity owing to the pres- 
ence of nodes or anisotropy, etc. While we can not find 
the decisive origin, the following discussion about the 
penetration depth suggests that the former possibility 
might be more likely in our film. o^iiS) can be rep- 
resented as 02(0;) ~ 1/ fiouj\ 2 (T) (the thick line in the 
upper panel of Fig. 2(b) is the lo^ 1 fitting curve), and 
we estimated the magnetic penetration depth at absolute 
zero (A(0 K) ~ 710 ± 70 nm) assuming the temperature 
dependence of the phcnomcnolo gical penetration depth 
is A(0 K) = \{T)yJl - (T/T c ) 4 . If we calculate the su- 
perconducting plasma frequency as uj ps = c/A, uj ps — 
0.28 eV. Because the carrier density is proportional to 
the square of the plasma frequency, this result suggests 
that 25% of the normal carrier contributes to the super- 
fluid density in our film. Therefore, we suggest that the 
former possibility about the finite conductivity in the su- 
perconducting state might be more likely, in the present 
stage. The penetration depth of our film is larger than 
the previously reported value, A(0 K) = 208 nm (from 
H c i—), A(1.7 K) = 245 nm (from muon spin rotation^), 
A(10 K) = 350 nm and A(0 K) = 360 nm (from optical 
measuremen t 23 ' 26 ) . While a slightly broad superconduct- 
ing transition (AT C ~ 2.9 K) might be also related to a 
larger A, there is a report for the classical superconduc- 
tor, lead, that A of films is approximately 10 times larger 
than that of bulk sample 27 . Therefore, we think that the 
interpretation for this issue deserve further studies, using 
various different films. We note that this will not affect 
the essential aspect of the following discussions regarding 
the superconducting gap. 

Next, we estimated the superconducting gap energy, 
A sc . At T < T c , the transmitted signal through the film 
is weak and noisy, so we extracted the data in the vicinity 
of main peak (4 ps < t < 11 ps) from the time-domain 
waveform. The upper panel of Fig. 3(a) shows the fre- 
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FIG. 3: (color online), (a) (upper panel) The fre- 
quency dependence of the transmittance of the ex- 
tracted time-domain waveform, (bottom panel) The 
normalized transmittance by the transmittance di- 
rectly above T c (20 K). The lines are the eyeguide, 
while the down arrows indicate the characteristic fre- 
quency where the transmittance increases in compar- 
ison to that of 20 K. (b) The temperature depen- 
dence of the superconducting gap energy, A sc , which 
is estimated from the characteristic frequency (down 
arrows in Fig. 3(a)). The solid line is the fitting 
curve of the BCS-type temperature dependence of 
A sc , A SC (T) = A tanh(1.74(r c /T- l) ' 5 ). 



quency dependence of the transmittance, T tr i m (cj), which 
is Fourier-transformed from the extracted time-domain 
waveform. By using this procedure, the noisy data in 
the original time-domain waveform (at t < 4 ps and t 
> 11 ps) does not affect the Fourier-transformed spec- 
trum such that we can discuss the physical nature in 
the superconducting state more clearly. Compared with 
Ttrim(w) just above T° nsct (T = 20 K), T trim (w) below 
j-onset wag ennancec j a t gome frequencies. The bottom 
panel of Fig. 3(a) indicates the normalized transmit- 
tance, T normalizod (iu,T) = T tlim (oj,T)/T tl i m (uj,T = 20 
K). The solid lines are the eyeguide. From this panel we 
can easily identify the characteristic frequency, lj* , where 
Tnormaiizcd^) = 1 (down arrows). Notably, uj* gradually 
increases with decreasing temperature, suggesting that 
hu* is related to the energy scale of the superconduct- 
ing gap. In dirty-limit BCS superconductors, Palmer 
and Tinkham calculated the frequency dependence of the 
transmission ratio2&. According to their calculation, Huj* 
is close to 1.4A SC , which leads to A SC (T = 0) — 2.0 meV 
in the case of our film. For strong-coupling superconduc- 
tors, however, Huj* ~ A sc . For the AeFe2As2 system, the 
strong coupling is suggested by a /xSR study-22. There- 
fore, we regard that fko* ~ A sc , and subsequently plot 
the temperature dependence of A sc in Fig. 3(b). The 
solid curve indicates the fit to the BCS-type function, 
A SC (T) = A sc (0)tanh(1.74(T c /T - l) 5 ), which is also 
applicable to strongly-coupled superconductors^. Using 
these fitting parameters (A sc (0) = 2.8 ± 0.5 meV and 
T c = 18 K) we obtained 2A sc (0)/fc B T c = 3.6 ± 0.6. We 
presume that the difference between T c determined by 
this fit and that obtained in the DC resistivity measure- 
ment originates from the slight change of quality in the 
sample handling in ambient air. According to the inset 
of Fig. 2(c), h/r is about 59 meV at T = 50 K. There- 
fore, we can conclude that our sample is close to the dirty 
limit (A sc (0) <C h/r). This result also supports the va- 
lidity of the above estimation of A sc for the dirty and 
strong-coupling material. Using ARPES measurement, 
Terashima et al} A reported that there are two supercon- 
ducting gap energies in a BaFei.85Coo.isAs2 single crystal 



(T c = 25.5 K). These gaps include the small supercon- 
ducting gap, Ag(0) ~ 4.5 meV, opens at the electron- 
type Fermi surface near the M point, and the large 
superconducting gap, A^(0) ~ 6.7 meV, opens at the 
hole- type Fermi surface near the T point. The value of 
2As(0)/fcBT c is close to the value we obtained, 3.6 ± 0.6 
(2A s (0)/fc B T c = 4.1 and 2A L (0)/k B T c = 6.1). There- 
fore, we concluded that the enhanced part of T tr i m (cj) at 
T < T c originates from the superconducting gap opened 
at the Fermi surface near the M point. The large su- 
perconducting gap opened at the Fermi surface near the 
r point exceeds our experimental frequency range (6.7 
meV corresponds to 1.63 THz). In completion of this 
work, quite recently an FIR measurement with using the 
Kramcrs-Kronig transformation^ reported <7i(w) indi- 
cating the existence of the superconducting gap at 3.1 
meV and 7 meV. Also, the optical measurement with us- 
ing the Mach-Zehnder arrangement^ reports the single 
small superconducting gap (1.85 meV). The supercon- 
ducting gap we obtained is consistent with that of opti- 
cal conductivity^, and hut* corresponds to the smaller 
superconducting gap. 

In conclusion, we measured the complex conductivity 
of Co-doped BaFe2As2 epitaxial thin film in the THz 
region, and we observed the Drude-type complex con- 
ductivity in the normal state. In the superconducting 
state, a(uj) of the film changes to the typical a(ui) of the 
superconducting material. We estimated the magnetic 
penetration depth at absolute zero as 710 nm. Also, we 
estimated the superconducting gap energy as 2.8 meV, 
which is considered comparable to the energy gap opened 
at the Fermi surface near the M point. We succeeded 
in obtaining the low-energy elementary excitation of Fe- 
based superconductor using the electromagnetic method 
without invoking Kramers-Kronig transformation. 
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